Mini-magnetospheres above the lunar surface and the formation of lunar swirls 
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In this paper we present in-situ satellite data, theory and laboratory validation that show how 
small scale collisionless shocks and mini-magnetospheres can form on the electron inertial scale 
length. The resulting retardation and deflection of the solar wind ions could be responsible for the 
unusual "lunar swirl" patterns seen on the surface of the Moon. 



Miniature magnetospheres have been found to exist 
above the lunar surface [1 and are closely related to fea- 
tures known as "lunar swirls" [2 j. Mini-magnetospheres 
exhibit features that are characteristic of normal plane- 
tary magnetospheres namely a collisionless shock. Here 
we show that it is the electric field associated with 
the small scale collisionless shock that is responsible 
for deflecting the incoming solar wind around the mini- 
magnetosphere. These ions impacting the lunar surface 
resulting in changes to the appearance of the albedo of 
the lunar "soil" [2] . The form of these swirl patterns 
therefore, must be dictated by the shapes of the colli- 
sionless shock. 

Collisionless shocks are a classic phenomena in plasma 
physics, ubiquitous in many space and astrophysical sce- 
narios [3 . Well known examples of collisionless shocks 
exist in the heliosphere, where the shock is formed by the 
solar wind interacting with a magnetised planet. What 
is a surprise is the size of the mini-magnetospheres, of 
the order of several 100 km; orders of magnitude smaller 
than the planetary versions. Results from various lunar 
survey missions have built up a good picture of these 
collisionless shocks. 

These collisionless shocks have a characteristic struc- 
ture in which the ions are reflected from a rather narrow 
layer, of the order of the electron skin depth c/u pe (where 
c is the speed of light and u pe is the electron plasma fre- 
quency), by an electrostatic field that is a consequence 
of the magnetised electrons and unmagnetised ions. The 
narrow discontinuity in the shock structure produces a 
specular reflected ion component with a velocity equal 
to or greater than the incoming solar wind velocity. The 
reflected ions from a counter-propagating component to 
the solar wind flow that form the magnetic foot region, 
which extends about an ion Larmor orbit upstream from 
the shock. This occurs when the Mach number (the ra- 
tio of flow velocity to Alfven velocity) is of the order 3 or 
less. 

We have carried out laboratory experiments using a 
plasma wind tunnel, to investigate mini-magnetospheres 




FIG. 1. The Reiner Gamma formation (7.4°N, 300. 9°E) is an 
example of a lunar swirl. Pictured here on the left hand side of 
the image. Reiner Gamma is named after the Reiner impact 
crater shown for comparison on the right. The crater is 117km 
to the east and has diameter of 30km with a depth of 2.6km. 
By contrast the unusual diffuse swirling of the formation and 
concentric oval shape has fluid-like "wisps" that extend fur- 
ther to the east and west. Its distinctive lighter colour stands 
out against the flat, dark surface of Mare Oceanus Procel- 
larum. Unlike crater ejecta, the tilted "7" shape of the for- 
mation appears unrelated to any topographic structures that 
would account for its presence. Image courtesy of NASA. 



and found that they show characteristics similar to the 
lunar mini-magnetospheres. A quantified comparison be- 
tween the observations, both in space and in the labora- 
tory, with theoretical values shows excellent agreement. 

The Reiner Gamma Formation shown in Figure [I] is 
one example of a number of small "swirls" of appar- 
ently lighter coloured material visible on the lunar sur- 
face. These distinctive patterns do not appear to corre- 
late with other surface features - such as impact craters 
or mountains and valleys, but do coincide with patches 
of significant magnetic field [2 . All of these lunar swirls 
have been found to be associated with magnetic anoma- 
lies [2] . The appearance of lighter albedo material on the 
Moon is usually indicative of the presence of younger or 
less weathered lunar material [4]. Explanations that ac- 
count for these changes in the albedo are either a) that 
the lighter regions have been "shielded" from receiving 
the same solar wind flux as the surrounding regolith and 
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FIG. 2. A sketch of the generic scenario of a miniature 
magnetic field emerging from the lunar surface and interacting 
with the solar wind. Flow is from left to right. 



hence appear younger [2 j, or b) that the lighter colour 
material originates from just below the surface and has 
been lifted up and deposited on the top of older/darker 
regolith [5]. The sharpness of the lunar swirl formations 
is enhanced by the contrast of "dark lanes" (suggesting 
locally enhanced solar wind proton bombardment) within 
the high-albedo swirls. 

Evidence by in-situ measurements from space craft 
including Lunar Prospector (1998-1999) [6 , Kaguya 
(2007-2009) 0, Chandrayaan-1 (2008-2009) [8 and No- 
zomi (1998) spacecraft [9 j, are consistent with the pres- 
ence of collisionless shocks and the formation of mini- 
magnetospheres. Because the observational data de- 
rives from a sequence of case studies from different 
missions [T[ ITU- 16 , it is to be expected that there is some 
variation in consistency. However a schematic picture of 
the interaction can be unravelled from the specific ob- 
servations to form a simplified, generic model. This is 
illustrated in Figure [2] 

A region of enhanced magnetic field strength (by fac- 
tors of 2 to 3) is observed at an altitude x si above the 
lunar surface where the magnetic intensity from the mag- 
netic anomaly reaches pressure balance with the plasma 
pressure from the solar wind. The magnetic field com- 
ponents were observed to rotate in a fashion consistent 
with the spacecraft passing through a region in which 
the solar magnetic field was being "draped" around a 
small magnetic obstacle or "bubble" pQ . Within the nar- 
row barrier region is a low density cavity seen in the ion 
data [T7J. The barrier region is of the order of kilome- 
ters across. Just ahead of a large ramp in magnetic field 
strength there is enhanced magnetic field turbulence in 
the solar wind magnetic field [1. This is accompanied 
by electrostatic solitary waves at the lower hybrid fre- 
quency and electron fluxes abruptly increase and their 
energy distribution changes indicating that electrons are 
energized and not simply compressed [1, 10-13, 15]. In- 
tense electrostatic waves, of frequencies of the order of the 
plasma lower-hybrid frequencies (l-10Hz) were recorded 



over the locations of the magnetic anomalies by the 
Kaguya spacecraft [12j [13] . Variations in the intensity 
of these waves over magnetic anomalies with solar wind 
pressure suggests a dynamic interaction [12]. At 100km 
Kaguya observed protons reflected back from the mag- 
netic structures with greater energies (by factors 3 to 6) 
than the incident solar wind flux [13]. Chandrayaan- 
1 also observed back streaming protons accelerated by 
similar factors close to the shock surface [TTJ [17] . These 
higher energy protons are accelerated by the convective 
electric field seen by the reflected protons in the solar 
wind flow [T3] , Below 50km altitude, proton back scat- 
tering disappeared [14], indicating that the flowing so- 
lar wind ions were no longer reaching these altitudes, 
and suggesting that the spacecraft was flying through 
a solar wind plasma "cavity". Chandrayaan-1 provided 
two-dimensional maps of the spatial extent of the cavity 
above the magnetic anomaly near the Gerasimovic crater 
[TT] . The dimensions of the cavity and the magnetic field 
anomaly (~360km across) were very similar. The cav- 
ity was more distinct in the higher energy ions >150 eV 
than at lower energies. The overall size was about 360 
km in diameter, coincident with the center of the mag- 
netic anomaly. An outer ring of about 300 km wide of 
enhanced proton flux suggested that the incoming solar 
wind ions were being deflected around from the central 
"bubble" to impact the surface in this relatively narrow 
surrounding region. Observations provided by the No- 
zomi spacecraft in 1998, suggest that characteristic wakes 
arising from the lunar mini-magnetospheres can extend 
to significant altitudes (2800 km). The instrumentation 
on-board recorded two peaks in density of non-thermal 
protons either side of a cavity of reduced ion flux [16]. 
These observation agree with those of Chandrayaan- 1 
at intermediate altitudes above the mini-magnetosphere 
bow shock [15] . 

The key to understanding how such structures can 
arise is to use a two fluid model of the plasma, in which 
the ions in the flowing plasma are un-magnetised and the 
electrons are magnetised. As the solar wind with its em- 
bedded magnetic field impacts a magnetic structure, a 
cavity is created. The cavity is bounded by an enhanced 
magnetic field that is about a factor three greater than 
the solar wind field with a width, L, estimated to be 
similar to the electron skin depth. This field traps a low 
density plasma that forms part of the barrier. The cavity 
is created by the induced currents that flow on its outer 
boundary, giving rise to magnetic field enhancement that 
opposes the penetration of the solar wind according to 
Lenz's law. The enhanced magnetic field drapes around 
the cavity. 

The magnetic field enhancement shown in the satellite 
observations, controls the flow of the electrons that are 
magnetized on these scales. The solar wind electrons are 
slowed and deflected by the magnetic structure. The ions, 
on the other hand, due to their greater inertia, cannot re- 
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spond sufficiently promptly to abrupt or sudden changes 
in magnetic field on these scales. These unmagnetized 
ions can therefore easily penetrate through the barrier. 
Ions flowing through the barrier result in a space charge 
separation between the electrons and ions, forming an 
electric field that is responsible for slowing and reflect- 
ing/deflecting of them over distances of the order of the 
electron skin depth [T8] . 

The equation of motion that controls the behaviour of 
the ions is 19 : 



d 1 2 



dpi 
dx 







(1) 



where is the electric potential ramp across the barrier 
and Pi is the ion pressure = riikTi and Vi is the ion ve- 
locity, the bulk value of which is the solar wind speed, 

The expression for the electric potential component, 0, 
responsible for slowing and deflecting the ions is 
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If the instantaneous density, n, here is ~ 5 x 10 6 m -3 and 
a magnetic field B z ~ 30 x 10 _9 T values similar to those 
observed at the pile-up reported by Lunar Prospector 
PQ, then the mean value from (2) of the potential would 
be (j)theory ~ 450V. This average value is very similar 
to the (pots ~ 400V [15] and would account for the ob- 
served counter- streaming protons. Although the force 
does act initially on the electrons, the resulting electric 
field formed then acts on the ions. The consequence is 
this force acts to keep particles out of regions of high mag- 
netic field. The shock thickness of the electric field, L, is 
much narrower than the ion Larmor radius. The ions ex- 
perience a sufficient impulsive force in a direction normal 
to the barrier to reverse their velocity. These ions move 
upstream with a velocity 2v sw in the solar wind frame 
and form a broad much thicker region than the barrier 
sometimes known as the shock foot. Within the shock 
foot region, the counter streaming ions are responsible 
for a number of micro-instabilities, such as the modified 
two stream instability [22 , that drives plasma wave tur- 
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bulence close to the lower hybrid frequency = (uj ce uj C i) 1 , 
where uo ce and uj C i are the electron and ion cyclotron fre- 
quencies respectfully. This would agree with data from 
all the spacecraft that observed intense lower hybrid elec- 
trostatic oscillations of the order of 1-10 Hz [16]. The re- 
flected or deflected ions can also form a non-thermal ring 
distribution in ion velocity space due to E x B pick-up 
exactly as reported [16] . Alternatively a simple counter- 
streaming population would be observed depending upon 
the particular conditions intersected by the spacecraft on 
its fly through. Both are consistent. The changes in the 
particle distributions observed by the in-situ spacecraft, 
between the ions and electrons either side of the shock re- 
sult from the formation and interaction of lower-hybrid 



TABLE I. A comparison of the absolute and dimensionless 
distances (normalized to c/ujpi), in space and the laboratory 
experiment. 





Space 


Laboratory 


Parameter in SI Units 


Value 


Dim' 


Value 


Dim' 


Thermal energy in eV 


5 




5 




Density, ra -3 


5.10 6 




10 17 




Flow speed, ms -1 


4.10 5 




8.10 4 




Magnetic field strength, T 


1 n — 8 
10 




0.03 




Plasma Beta 


0.1 




2.10" 4 




MACH# Acoustic, Alfvenic 


20, 5 




3.5, 0.03 




Flow Mean free path, m 


10 16 


10 11 


300 


400 


Debye length, m 


7 


10" 5 


5.10" 5 


7.10" 5 


Electron Larmor radius, m 


800 


8.10" 3 


3.10" 4 


4.10- 4 


Thermal Ion Larmor radius 


6.10 4 


0.6 


10" 2 


0.01 


Flow Ion Larmor radius, m 


5.10 5 


5 


3.10" 2 


0.04 


Electron Skin depth, m 


2.10 3 


0.02 


2.10" 3 


0.01 


Ion Skin depth c/ujpi, m 


10 5 


1 


0.7 


1 



waves generated close to the bow shock [23]. The res- 
onant interaction between lower-hybrid turbulence and 
electrons can result in field aligned electron acceleration 
[23] . These waves are most probably excited by the mod- 
ified two stream instability driven by reflected ions. 

The density of transformed electrons, n^, is estimated 
by balancing the growth rate of the instability initiated 
by pickup ions with Landau damping due to electrons 
moving parallel to the magnetic field. Estimations for 
the average energy e e of the accelerated electrons and 
their number density n^ e can be made [23]: 



v 2 / 5 



v 2 / 5 



1/5 
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(3) 



(4) 



where a is the energy transformation efficiency from ions 
to electrons (with masses mi and m e respectfully). 

The energy transformation coefficient a, has only a 
weak influence on the result [23] . Therefore using a value 
of a ~ 0.1, the case near a bow shock, with an ion en- 
ergy of ~ IkeV and electron energy e e ~ lOOeV, the den- 
sity of the accelerated electrons would be 10% of the ion 
density (n^e — O.ln^). This is consistent with the typi- 
cal values reported by Kaguya [14] and Chandrayaan-1 
[24] . This simple estimate demonstrates the large effi- 
ciency of lower-hybrid waves as an acceleration mech- 
anism for the electrons. The result is a more efficient 
boundary than would be predicted by MHD (magneto- 
hydrodynamics); comparisons with the particle distribu- 
tion data could confirm this. The theoretical width of 
the barrier is expected to be of the order of the electron 
skin depth, L =1 to 2 km, and not the 100's km of the 
ion skin depth [18 . See Table [i] for typical values. 
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FIG. 3. Photographs, with graphical labels overlaid, of the 
supersonic plasma stream being deflected by two different 
strength magnets, (a) 23x20mm, 0.45 T on axis. (b)10x3mm, 
0.30T on axis. External horizontal magnetic field is 0.03T. A 
field-aligned current can be seen connecting the 'cusp' region 
to the southern poles in both cases. 
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FIG. 4. In-situ Langmuir probe plots of the plasma den- 
sity and floating potential ahead of the larger target magnet 
(23x20mm, 0.45T). Data from [26]. 

Computational simulations of kinetic processes are 
non-trivial. A scaled laboratory experiment has the po- 
tential to deliver physical insights and to offer observa- 
tional signatures which can be examined, to establish 
whether they are consistent or not with the in-situ space 
data. The Plasma Wind Tunnel [25] used here shares the 
same phenomenological regime as the subset of analogous 
space parameters. These are: the plasma is collisionless, 
the bulk flow speed is supersonic and the electrons are 
magnetized but the ions are not. The absolute and di- 
mensionless parameter are given in Table [IJ 

In Figure [3] we show two photographs of the Plasma 
Wind Tunnel in action. The supersonic hydrogen plasma 
stream is here impacting the magnetic dipole fields of 
two different sized magnets under the same conditions. 
The plasma stream appears dark on the photograph and 
is enhanced for clarity. The structures are 3D so some 
blurring due to line-of-sight occurs. Despite the > x2 
difference in the size of the magnetic "bubble" obstacle 
in the two cases, the plasma stream is corralled the same 




FIG. 5. A plot of the measured ion density (black dots) and 
electric potential (green triangles) along the axis of the ex- 
periment in the vertical plane. The plasma flow is from left 
to right. Data from [26 j. 

way into a narrow boundary layer and deflected around 
the outsides to form a cavity. The key inset shows the 
ion Larmor radius to scale with the photographic data 
demonstrating that the interaction is on the sub-ion Lar- 
mor radius scale as with the Lunar case. 

The Langmuir probe data shown in Figure [4] corre- 
sponds to the larger of the two magnets shown in Fig- 
ure [3J The ion density shown in the upper panel con- 
firms the formation of a cavity or void in the plasma 
stream. The floating electric potential (lower panel) co- 
incides with the locations of the measured ion density 
shown in the upper panel. This demonstrates that the 
ions are indeed electrostatically confined on this scale. 
The measured width of the boundary layer is between 2 
to 3mm, similar to the calculated electron skin depth of 
~2mm [18] and is very much less than the calculated ion 
inertial length of 700mm or the 30mm of the ion Larmor 
radis (see Table [l]). The quantified relation between the 
ion density and the electric potential can be better com- 
pared in the ID plot shown in Figure [5j At the point 
x s shown in the figure, the measured electric potential is 
$obs ~ 8V. Here again, the estimate provided by Equa- 
tion (2) provides a good match to the observed value with 
^theory ~ 12.5V (using B = 0.05T, n = 5 x 10 16 m- 3 ). 
This shows that the dimensionless analysis holds on a 
very different absolute scale to the case in space. 

The analysis presented here shows that the strength 
of the deflecting electric field of a mini-magnetosphere 
and collisionless shock, is not dependent on the overall 
size of the magnetic "bubble" but is related to the local 
gradient in the magnetic field strength. None of the real 
features are simple single dipoles. Close to the surface, 
the magnetic topology in a magnetic anomaly is likely 
to be very irregular, comprising a range of overlapping 
cavities and gradients. This would lead to a pattern of 
retarded and accelerated space weathering and hence ar- 
eas of lighter material with embedded "dark lanes". A 
further range of contrasts in the surface "aging" would 
come from the variations in the impacting plasma wind 
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environment with solar activity, the Moon's orbit in and 
out of the Earth's magnetosphere and lunar phases. 

In conclusion, the model of small scale collisionless 
shocks that we present agrees with a laboratory plasma 
wind tunnel experiment confirming the presence of a nar- 
row electrostatic potential of thickness of the order of the 
local electron skin depth, and that it is the force respon- 
sible for the control and deflection of the ions. All the 
observational data from spacecraft is quantinably consis- 
tent with the theoretical model. 
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